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By combining effective field theory with the standard nuclear physics approach (SNPA) we obtain a 
high-precision estimate of the S factor for the solar hep process. The accurate wave functions available 
in SNPA are used to evaluate the nuclear matrix elements for the transition operators that result 
from chiral perturbation theory (ChPT). All the contributions up to N 3 LO in ChPT are included. The 
resulting parameter-free, error-controlled prediction is: S(hep) = (8.6 ± 1.3) x 10~ 20 keV-b. 

This brief report is based on the results of work done in collaboration with L.E. Marcucci, R. 
Schiavilla, M. Viviani, A. Kievsky and S. Rosati || |3j. A detailed exposition of the basic ideas 
! underlying our approach can be found in Jl|] . 
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O . 1 Introduction 
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The hep process in the Sun 



O ' 3 He + p -> 4 He + e+ + v . (1) 

produces the highest energy solar neutrinos, S™ ax (/iep)~20 MeV. The hep neutrinos therefore can 
influence the interpretation of the results of a recent Super-Kamiokande experiment that have 
raised many important issues concerning the solar neutrino problem and neutrino oscillations J|,[5| . 
It is to be noted that the reliable estimation of the hep cross section, indispensable for addressing 
these issues, is a long-standing challenge for nuclear physics §. This is mainly because the leading 
^ ■ one-body contributions are highly suppressed and furthermore the chiral filter mechanism - which 

allows us to accurately estimate many-body corrections - is ineffective for this process. For a detailed 
discussion, see Ref. [@,||. 

The objective of our present work is to obtain a significantly improved estimate of the hep 
rate using effective field theory (EFT). To this end, we adopt a strategy that exploits the known 
merits of the standard nuclear physics approach (SNPA) and heavy-baryon chiral perturbation 
theory (HBChPT) simultaneously. HBChPT, a well established low-energy EFT of QCD, is used 
to calculate the transition operators; all the operators up to next-to-next-to-next-to-leading order 
(N 3 LO) will be considered. The evaluation of the corresponding nuclear matrix elements requires 
highly accurate nuclear wave functions. Although it is, at least in principle, possible to derive from 
HBChPT nuclear wave functions to a specified chiral order, we choose not to do so. Instead, we 
use realistic wave functions obtained in SNPA. The power of the proposed scheme is the ability 
to correlate the beta decay processes in the A = 2,3,4 nuclei. Thus, as explained in more detail 
below, if the single unknown constant in our EFT is fixed using one of these processes, then we can 
make totally parameter-free predictions for the remaining processes. 
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2 Theory and Results 

The hep process is dominated by the Gamow- Teller (GT) transition, and hence the reliability of 
the hep rate estimate is essentially governed by precision with which one can calculate the GT 
amplitude. According to the chiral counting rule Q, the leading order contributions are due to 
the well-known one-body (IB) operators, and the first corrections arise from N 3 LO two-body (2B) 
currents that are suppressed by (Q/A x ) s compared to the IB. Here Q stands for the typical three- 
momentum scale and/or the pion mass, and A x ~ 1 GeV is the chiral scale. As stated, we consider 
here a^the contributions up to N 3 LO. It is worth emphasizing that three-body currents, which are 
N 4 LO, can be legitimately ignored in our N 3 LO calculation. 

The N 3 LO 2B currents consist of the one-pion-exchange (OPE) and nucleon-nucleon contact- 
term (CT) parts, A2B = A2b(OPE) + A2b(CT). With the low-energy constants fixed from irN 
data |i~3|l ) the OPE part is completely determined. On the other hand, the CT part contains one 
parameter, d R , whose direct evaluation from QCD is not available at present. Fortunately, it turns 
out that tritium /3-decay, //-d capture and v-d scattering are sensitive to the same parameter, d R , 
and that they do not depend on any additional parameters up to N 3 LO. Thus, any of these processes 
can give the renormalization condition to fix the value of d R . Here we choose to use the tritium 
/3-decay rate, Tp, which is accurately known experimentally [jl4). Once d R is fixed, our calculation 
involves no unknown parameters. 

We calculate the matrix elements of the transition operators with state-of-the-art realistic 
nuclear wave functions. We employ the correlated- hyperspherical-harmonics (CHH) wave functions, 
obtained with the Argonne v\% (Avl8) potential (supplemented with the Urbana-IX three-nucleon 
potential for the A > 3 nuclei) 0. To control short-range physics in a consistent manner, we apply 
the regulator 

Sa(q 2 ) = exp f- Jtlj . (2) 

to all the nuclear systems in question. The cutoff parameter A characterizes the energy-momentum 
scale of our EFT. 

The value of d R determined from the experimental value of Tp is d R = (1.00 ± 0.07, 1.78 ± 
0.08, 3.90 ±0.10). Here and hereafter, parenthesized three numbers correspond to the three choices 
of A, A = 500, 600 and 800 MeV, in this order. To see the role of the tr^-term, it is informative 
to look at <52B> the ratio of the 2B matrix element to that of IB. With only the OPE part taken 
into account, we have 5^ E = (-1.1, -1.5, -2.0). The inclusion of the d R term leads to 5^ LO = 
^m E+CT = ( — 0.60, —0.64, —0.73). Thus the J R -term drastically reduces the A-dependence of 
the 2B contribution; we see only ~10 % variation for the entire range of A under study. The 
A-dependence in the total GT amplitude becomes more pronounced due to a strong cancellation 
between the IB and 2B terms, but this amplified A-dependence still remains within acceptable 
levels. 

In addition to the 3 Si contributions governed by the GT amplitude, there are also tiny 1 S'o 
and sizable P-wave contributions. The latter have little A-dependence (< 2 %), and responsible for 
about one-third of the total S factor. Adding all the contributions, the S-factor at threshold reads 

S(hep) = (8.6 ± 1.3) x 10~ 20 keV-b , (3) 

where the "error" spans the range of the A-dependence for A = 500-800 MeV. This result is to be 
compared to the latest SNPA estimate in Ref. @: S = 9.64 x 10~ 20 keV-b. 
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3 Discussion 

By determining the only parameter of the theory, d R , from the experimental data on triton beta 
decay, we have succeeded in estimating the hep S-factor in a parameter-free and error-controlled 
manner. Our HBChPT calculation (up to N 3 LO) gives a much more accurate estimate than hitherto 
available. The EFT results turn out to give support to those obtained in the latest SNPA calculation 

0. 

To decrease the uncertainty in Eq.(^), we need to reduce the A-dependence in the two-body 
GT term. According to a general tenet of EFT, the A-dependence should diminish when higher 
order terms are included. A preliminary study indicates that it is indeed possible to reduce the 
A-dependence significantly by including N 4 LO corrections. In this connection, we remark that the 
hen process, 3 He + n — > 4 He + 7, seems very interesting to look into. The hen process shares many 
features with the hep including the suppression of the IB matrix element and the structure of the 
many-body currents. Accurate experimental data are available for the hen cross sections, but so 
far no theoretical calculations have succeeded in explaining the data quantitatively. Thus applying 
the same EFT technique to the hen process |l5| is very interesting, and that will also provide a 
useful check of the formalism employed in our estimation of S(hep). 
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